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ABSTRACT 


Electromagnetic pulse has become a possible threat to 
nearly all sophisticated military systems. The crossed di- 
pole receiving antenna has been used as a representative 
model to approximate electromagnetic pulse effects on air- 
craft. A very basic way to approximate an aircraft structure 
is by a two-dimensional cross. 

This paper iS an application of the Piecewise-Sinusoidal 
Reaction Matching Technique (PSRMT) to find the current den- 
sity distribution in a two-dimensional cross illuminated by 
a monochromatic plane wave. Other two-dimensional structures 
are solved previous to the cross structure in order to gain 
insight and to validate this approach with respect to pre- 


vious solution techniques. 
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i INTRODUCTION 


pee WEED FOR STUDY 

Electromagnetic pulse has become a possible threat to 
nearly all sophisticated military systems [Ref. 1]. Under 
the proper circumstances a significant portion of the energy 
released during a nuclear detonation will be appear as an 
Electromagnetic Pulse (EMP) having the same frequencies as 
those employed by most of our commercial radio and military 
system equipments. 

Any metallic object, such as an aircraft structure, ex- 
posed to electromagnetic fields can be a collector of elec- 
tromagnetic energy, that is, act like an antenna, even 
though it was never intended to be that. 

The crossed-dipole receiving antenna has been used as 
Emerprcsencative model to approximate electromagnetic pulse 
effects on aircraft [Ref. 2]. The crossed-dipole, in this 
case the aircraft, can be modeled in its simplest way as a 
two-dimensional structure, that is, a cross with thickness 


equal to zero. 


B. BACKGROUND INFORMATION 

Electromagnetic boundary-value problems can be solved 
exactly via classical separation of variable analysis only 
for a few geometries, such as the sphere, spheroid, circular 
cylander, elliptical cylinder, strip and wedge. Such solu- 
tions can be expressed in terms of a summaticn of a set of 


eigenfunctions which can be evaluated with a high speed 


ale 





computer. However, for problems involving complicated geo- 
metries the exact solution is not available. Thus, the ap- 
proximate and numerical methods are of great significance. 
The aeseerech discussed herein is directed toward the numeri- 
cal solution of the problem of scattering from conducting 
bodies of arbitrary shape. 

Two methods [Ref. 3] are available for electromagnetic 
modeling of continuous conducting surfaces with arbitrary 
shape: the wire grid model and the surface current model 
uSing rectangular-pulse bases. Both methods have similar 
limitations with the maximum cell width restricted to ap- 
proximately A/10. Unless the conducting body is symmetric 
or is a figure of revolution, computer storage requirements 
have limited the moment-method application to bodies with | 
surface area not exceeding one or two square wavelengths. 

In this paper, the piecewise-sinusoidal reaction tech- 
nique applied by Richmond [Ref. 3] to scattering by conduct- 


ing bodies of arbitrary shape is used. 


Seo lo EMENT OF THE PROBLEM 

The objective of this research is to find the current 
distribution in a two-dimensional cross. In order to get 
an understanding of the "junction problem" that will be 
present in the cross structure, the induced current density 
on a perfectly-conducting rectangular plate will be first 
solved for, this result having been presented by Richmond 


and Wang [Ref. 3]. Likewise the square plate solved by 


Rahmat-Samii and Mitra [Ref. 4] will be studied. From this 
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the complexity of the structure will be increased, solving 
for current distribution on structures like an L, T and 


finally the cross. 


Bees SCOPE AND LIMITATIONS OF THE STUDY 

The general solutions steps are: 

1) Divide the continuous conducting surface into cells 
with the surface current distribution being expanded in 
overlapping sinusoidal bases functions. 

2) Via an application of Galerkin's method, the inte- 
gral equation formulated with the zero-reaction concept is 
reduced to a matrix equation. 

3) The current distribution over the conducting surface 
1s determined via matrix inversion. 

The manner in which the structure is subsectioned is 
very important in any numerical analysis approach. Likewise 
the basic current distribution chosen for the subsections 
drastically affects the amount and difficulty of calculations. 
These methods will yield approximate results, the quality of 
the approximation depending on the amount of effort (and 
computer time) one is willing to devote to the problem. 

As will be seen symmetry may simplify considerably the 
number of calculations involved. In some cases one may want 
to take advantage of symmetry; in other cases it may be 
Simpler to ignore it. It appears that one must trade some 
formulation time for computer time. Is the tradeoff worth- 
while? The answer to this question will depend on the extent 


of the computations. In reducing the number of unknowns, one 


Lu 


increases the analysis time but reduces computer programming 


time and computer running time. 


Ls 





Pi. METHOD OF ANALYSIS 


Pee laOD OF MOMENTS 
The method of moments has been widely explained in the 
literature. The interested reader may refer on this subject 


to References 5, 6 and 7. 


EyeeePLECEWISE SINUSOIDAL BASIS SETS 

The physical problem may be represented by many different 
operator equations, and a suitable one must be chosen. There 
are also an infinite number of sets of expansion functions 
and testing functions that may be chosen. 

Intuitively, the more nearly the expansion set resembles 
the actual form of the current distribution, the fewer the 
terms required to accurately represent the current. In ad- 
dition, certain types of basis functions may reduce many 
time-consuming numerical operations incorporated in approxi- 
mations. The piecewise-sinusoidal basis set exhibits these 
qualities. The objective is to represent the current dis- 
tribution, over individual subsections of a domain with ex- 
pansion functions that possess greater pense fe Cine potential 
than do flat pulses and yet lead to well-conditioned matrices 
in approximate solution procedures. The application of piece- 
wise-sinusoids in Galerkin's method and the methed of moments 
to electromagnetic field problems is known as the Piecewise- 


Sinusoidal Reaction Matching Technique. 





C. PIECEWISE-SINUSOIDAL REACTION MATCHING TECHNIQUE. SURFACE MODEL 

The reaction concept and its application has been dis- 
cussed by Rumsey [Ref. 8], Cohen [Ref. 9], Harrington [Ref. 
10] and Richmond [Ref. 11]. 

The following material has been previously presented 
by Richmond and Wang [Ref. 3]; it is presented here for 
completeness. 

Consider the exterior scattering problem illustrated 
in Figure 1. In the presence of a dielectric or conducting 
body, the impressed electric and magnetic currents (J; 5M.) 
mopeeore the electric and magnetic field intensities (E,H). 
For simplicity let the exterior medium by free space. 

Eerom the Pin raeeeeaut valence EicomemlOn Semel kun@rs 5 
fiiementerior field will vanish (without disturbing the ex- 
terior field) if the following surface-current densities 


eee introduced: 


J. =n xH (1) 
ue =e ect (2) 
on the closed surface s of the scatterer. (The unit vector 


n is directed outward on s.) In this staon. illustrated 

in Figure 2, it is possible to replace the scatterer with 

free space without disturbing the field anywhere. 
Byedetinition, the incident field (E; ,H.) is generated 


by (J.,M.) in free space, and the scattered field is 


EE. ee E. (3) 
Eee) Her (4) 


diy 





J. 4M. 
—) — 
SOURCE 





oCATTERER 


Le 


Figure 1. The Source (J. M. pecenet@atcomene Pteld (Eon )s with 
scatterer. 





Js Me 
SOURCE 


(5) 
See ERER OK 
Plea ls Seyler 


(E,H) 
(JM 
Figure 2. The Interior Field Vanishes When the Currents 
(J_,;M_) are Introduced on the Surface of the 


See ee ne. 
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When the surface current (J. »M_) radiates in free space, it 
generates the fieid (E.,H.) in the exterior and (-E. ,-H.) 
in the interior region. This result, illustrated in Figure 
3, 1s deduced from Figure 2 and the superposition theorem. 
With the scatterer replaced by free space, it can be 
noted in Figure 2 that the interior region has a Nnivide fae ie 
As shown in Figure 4, by placing an electric test source J. 
in this region it is found from the reciprocity theorem 
that 


f (J°E,-M.cH,)ds + SSS(d,°E,-M,*H,)dv = 0 (5) 


where (EH) is the free space field of the test source. 
In words, Equation (5) states that the interior test source 
has zero reaction with the other sources. This "zero reac- 
tion theorem" was developed by Rumsey LRef. 8]. 

Equation (5) is the integral equation for the scattering 
problem, and the objective is to use this equation to determine 
the surface current distribution J. and M.. To accomplish 
this, these functions are expanded in finite series yielding 
a finite number N of unknown expansion constants. Next N 
simultaneous linear equations are joer to permit a solu- 
tion for these constants. One such equation is obtained from 
Equation (5) each time a new test source is set up. 

If the scatterer is a perfect conductor, the magnetic 
current M. vanishes. Assuming a finite conductivity and 


using the impedance boundary condition 


x2 el (6) 


le 





—_ = = e288 «= oe = 
—_ = ° — a 
—_ 


5 (E. Ho) 


Figure 3. The Exterior Scattered Field may be Generated by 
(J. »M_) in Free Space. 
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Meso eenmehieectrie Lest Source J, is Positioned in the 
Interior of the Scattering Region. 
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where Ze denotes the surface impedance equation (5) yields 
- e ee oa vd.e -_- ~? 
# J. LE, (nxH, )Z_ Ids fFICT. EB. M, +H, )dv. CD 


For two-dimensional problems involving arbitrary scat- 
terers, J. and M. are functions only of the position on the 
surface of the scatterer. If M,. vanishes and denoting 
(EH) as the free-space field of test-source m, Equation 


(7) yields 


piel GE (nx )2Z Ids = ff J, -E ds. (8) 


Now representing the electric current distribution as 
follows: 


N 
See a) (9) 


Js ln -—-n 


where the complex constants t, are samples of the function 
J,- The vector functions J. are known as basis functions. 
Basis functions J: and testing (weighting) sources J with 
unit current density at their terminals are used. 


From Equations (8) and (9) the following simultaneous 


linear equation is obtained 


N 
Tamia e mee =! hk : (10) 


teem = 1,2,3,...,N, where 


cae ak Ue ees (ibd) 2 d= 


= -Sf, Ve Oise ds (11) 
fie = ff. JE ds = ff JE: ds (129 


Zi 





The integrations in Equations (11) and (12) extend over 
the region where the integrand is non-zero. Region "m'" 
covers the interior test source, J? and region "n" is that 
portion of the surface s covered by the expansion function 
Jo: The first and the second integrals in Equation (11) are 
related by the reciprocity theorem, where EY is the free- 
space field generated by J. and the associated magnetic 


current Me 


D2 UNIDIRECTIONAL CURRENT MODEL 

A two-dimensional surface dipole located on the y-z plane 
1s illustrated in Figure 5(a). This source is an electric 
surface-current density with height 2D and width w. The 


surface current density is given by 


J) = sank One for 0 < 7 < 8 (13) 
- sink(Dtz) ~ Gorm= 05 <— zo 0 (14) 


Te) 7 sinkKD 


where K = 21/d. 
As illustrated in Figure 5(b), the current density is as- 
sumed to be piecewise-sinusoidal in the z-direction, with 
period equal to A, vanishing at the eee zZ = Dand z= -D, 
and equal to unity at the center. From Figure 5(c) it can 
be seen that the current density is assumed to be uniformly 
distributed in the tranverse direction. 

iiemempoterctonent distribution, Equations (13) and 
(14), will be used as the basis functions (J. ifipe oe] We 1 Ona 
(2)) for expanding the unknown current distribution induced 


on a conducting surface. Also, surface dipoles will be 


a 





O>0,D) (050,53) 
(00,0) . 
9V 5 y (0 ,0,0) Siesy 
Z 
oo,0,-D) (0,0,-D 
(a) (b) 
ay) 
uc 
0 W y 
Ce) 


Figure 5. An Electric Surface Dipole and Its Current 
Ptoiiay Distr bution . 
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employed as test sources with the reaction concept to solve 


the integral equation. 


Peecalewlattonsor the Coupling Term 


A 


The complex number C om in Equation (11) represents 
the mutual coupling between two sources. Let <m,n> denote 
the mutual coupling between sources, then Equation (11) can 


be written as 


Can = =< = ~ff J EL ds ClES3 


Relating the mutual coupling between two sources to the cir- 


cuit parameters [Ref. 3] 


C = V iL Cre) 


where ae 1s the open circuit terminal voltage induced at 
source m by source n, and Ton 1s the short circuit terminal 
current induced at source n. 

The surface dipole self impedance has been defined [Ref. 


Saleh the induced-emf formulation: 


_ mm _ mm 
Ls 7 ee C17) 
om om om 


Z a Hs eles 0) 2 ds. Cis) 


where J 1s the surface-current density of source m and En 
is its free space electric field. 
The mutual impedance between two surface dipoles is 


defined [Ref. 3] by 


2u 


So 
Zo = eo SL 


ae, TT we aS. (19) 
Gmnon 


J 
1 aio) 


The evaluation of the integrals in Equations (18) 
and (19) requires use of numerical integration techniques. 
K. Khemayodhin has developed a computer program to evaluate 
these integrals. Details and results are shown in his 
Engineer's Degree Thesis to be presented to the Electrical 
Engineering Department, Naval Postgraduate School, in 
December 1974. His computer program has been used in this 
thesis to calculate the self and mutual coupling between 
surface dipoles. 

Table I lists the self impedance of a center-fed 
two-dimensional sliipaile as a function of size. Table II 
lists the mutual impedance between the two-dimensional di- 
poles shown in Figure 7 as a function of their relative 
Besttion. 

Hrenmeequations (l/)Jeand (19) the mutual coupling 


memaerived as 


Cee ee Zo | (20) 


For the dipole shown in Figure 5, I is the total 


om 

Current at z = 0. Since J = 1 at the center of the dipole 
W 

an J dy = w. (21) 


For identical dipoles and from Equation (21), Equa- 


men (20) yields 


€ = ee re (oD) 





Table I. Self Impedance of Center-Fed Two-Dimensional Dipole 
Shown in Figure 6. 






a 


“ee 


Figure 6. Two-Dimensional Dipole. 
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Table II. Mutual Impedance of Center-Fed Two-Dimensional 
Dipoles Shown in Figure 7. 
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ee 


Figure 7. Coupled Two-Dimensional Dupotccewatin » = 0.25) 
and w = 0.25). 
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Came mieulatron of the Excitation Column 
In the matrix equation ae I = A? the complex 
quantities A. form the excitation column. Physically, co 
1s the reaction between the impressed source and the surface 


dipole m. If the impressed source is an electric line 


Seuree with current density J. Equation (12) reduces to 


°E. ds. (23) 


An = Tin Sn" Bi 


m. 


Assuming a line source located at a great distance from the 
surface, the incident field may be regarded as a plane wave 
with 


a Elles (2) 


For Eo = 1 and at x = Q, Ee = 1. For the two-dimensional 


dipole shown in Figure 5(a). Equation (23) reduces to 


W al 
AL f J dz dy. (25) 


O Oo 


Using the surface current density assumed in Equation (13) 


and (14), Equation (25) yields 


_ 2 W 
An Bes a ICD (1 - cosKD). (26) 


Peo GoNERAL CURRENT MODEL 

A two-dimensional surface quadrupole located on the y-z 
plane is illustrated in Figure 8(a). This source is an elec- 
tric surface-current density with height 2D and width 2 w. 


The surface current density is given by 
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. sin K(w-y) 
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where K = 21T/X. 

As illustrated in Figure 8(b) and (c), the current den- 
Sity 1s assumed to be piecewise-sinusoidal in the y-direction 
and z-direction, with period equal to A, vanishing at the 
edges and equal to unity at the center. Comparing this struc- 
ture with the one used in the Unidirectional Current Model 
1t can be appreciated that the quadrupole is a combination 
of the dipoles used in the former. Its current distribution 
will be the product of the two functions shown in Figures 
8(¢b) and 8(c). 

K. Khemaydhin has also developed a computer program to 
evaluate the current distribution on conducting bodies of 
arbitrary shape. Details are contained in his aforementioned 
thesis. This computer program has been used in this thesis 
to calculate the current distribution in two-dimensional 


structures when current in two directions was assumed. 
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ie COMPULER PROGRAM. DISCUSSION, DESCRIPTION AND RESULTS 


This section is concerned with the discussion and de- 
scription of the computer programs used to find the current 
distribution in the several two-dimensional structures to 
be treated. There are a number of refinements which can be 
introduced to improve the accuracy of the sclution, as well 
as to increase the effectiveness of the use of the computer 
memory. The principles involved are simple in concept; in 
each case an integral equation is replaced with a matrix 
equation, which is inverted to solve for the unknown source 


mime: On. 


A. UNIDIRECTIONAL CURRENT MODEL 

In the solution for current distribution in the several 
two-dimensional structures in which only one direction of 
current was allowed, it was assumed that a line source is 
located at a great distance from the surface so that the 
incident field may be regarded as a plane wave. The plane 
Peljeets assumed to be incident at broadside and in the 
z-direction. 

The number of subsections is, of course, a very important 
quantity in any Numerical Analysis Approach. One of the 
easiest methods of improving accuracy iS merely to increase 
the number of subsections, but increasing the number of sub- 
sections increases the number of unknowns, the time of com- 
puter programming, the computer running time and the computer 


memory space. Thus it appears that it is a tradeoff between 
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the extent of computations and the accuracy of the solution 
desired. 

For all the structures solved, the surface was divided 
into an integer number of equal two-dimensional dipoles. 
Although this method of segmentation gives an equal resolu- 
tion at the sharp edges (where the current density is vary- 
ing very rapidly) as well as in the center of the surface 
(where the current density is almost constant), it allows 
the use of symmetry for the calculations of the coupling 
terms. The use of symmetry, as explained in each of the 
particular cases, has saved considerable computational time. 

Mee liceekhectangulan Plate 

The first problem solved was the rectangular plate, 
for which Richmond and Wang [Ref. 3] have previous results, 
The rectangular plate, as shown in Figure 9, was divided 
into ninety overlapping dipoles. The current density dis- 
tribution was found with Computer Program 1 listed in Ap- 
pendix A. Its description is as follows: the program 
written for a rectangular shaped structure is a general pro- 
gram applicable for different methods of segmentation. The 
only restriction is that the segments Ginoiee) must be 
equal. The number of dipoles is limited by the size of the 
matrix which the inverting matrix subroutine can handle, 
usually 100. Subroutine CMTRIN, modified to handle matrices 
up to 200 by 200, was used in this program. 

The way to describe the rectangular plate is by en- 
tering the following parameters: (a) Ll is the number of 
Gapoles in a row (30 in Figure 9); (b) L3 is the number of 
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Hows Of dipoles (3 in Figure 9); (c) L7 is the total number 
of dipoles in the surface (90 in Figure 9); (d) D is the 
half-length of the dipole (see Figure 5(a)); (e) w is the 
width of the dipole (see Figure 5(a)). 

DO LOOPS 100, 200, 300, and 400 generate the coor- 
dinates of the source dipole and the test dipole as shown 
in Figure 7. Each dipole is at one time a test dipole, 
with all the other dipoles becoming source dipoles. In 
calculating the self impedance, the Complex Function ZSDPAR 
provides a 0.001A separation between dipoles and calculates 
this closely-coupled case of mutual impedance as the self 
impedance. 

The number of times the calculation of the coupling 
between dipoles has to be performed depends on the number 
of dipoles on the surface. The time each calculation of 
the coupling term takes in the IBM 360/67 computer has been 
estimated to be an average of 1.6 seconds. To perform the 
Seumealculations for the rectangular plate would take 216 
minutes of computer time. However, noticing that the coup- 
ling between dipoles depends only on their relative posi- 
mem, it is possible to Epp ly some symmetry to the problem. 
Pamewlating only the first row of the coupling matrix, that 
1s, the coupling between dipole 1 and all other dipoles in 
the surface will yield sufficient terms for filling the rest 
@eeene Matrix. So, only 90 calculations of the coupling 
term need to be performed, reducing the computer time dras- 


tically to 2.4 minutes. DO LOOPS 450 and 466 perform the 


34 





symmetry comparison by comparing the relative position of 
the dipoles. 

The calculation of the right hand side vector of 
Equation (10) is performed with the closed form of Equation 
C26). 

The coefficient matrix is then inverted using sub- 
routine CMTRIN. 

DO LOOPS 2400 and 2500 perform the multiplication of 
the inverted matrix [C(m,n)]7+ by the excitation column 
[A(m)]. 

DO LOOP 3000 performs the calculation of the current 
density at each dipole and normalizes it to the Magnetic 


Field Intensity 


Be 
n= aa C2) 
y 
From Equation (24), E, = 1, then 
ab 
Fe = = Se 
y n — 


The absolute value of the current density is plotted 
using the three-dimensional plot sub-program PLT3D1 and is 
shown in Figure 10. 

wee ie Sauare Plate 

The two-dimensional square plate is shown in Figure 
li. The structure was divided in 90 overlapping dipoles. 

The current density distribution was found using 
Computer Program 1 following exactly the same procedure as 
in the rectangular plate. The only difference is in the 


size of the dipoles. For the square plate D = 0.25,A and 
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Rectangular Plate (1A by 0.5) Magnitude of 


Induced Current. 


Figure 10. 
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Figure 11. Square Plate Segmented in 90 Overlapping Di- 
poles. Plane Wave Incident at Broadside. 
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w = 1/30A, instead of D = 0.1254 and w = 1/30A as used in 
the rectangular plate. 

The current density distribution is shown in Figure 
12 and was plotted using the three-dimensional plot sub- 
Program PLT3D1. 

3. ‘L'=Shaped Structure 

The 'L'-shaped structure is shown in Figure 13. The 
structure was divided into 60 dipoles. 

Computer Program 2 listed in Appendix A was used to 
solve the current density distribution for the 'L'-shaped 
structure. This is a general program applicable for dif- 
ferent ways of segmentation of an 'L'-shaped structure. It 
has the same restrictions as Computer Program 1 (rectangular 
piate), in fact both programs are quite similar. Combvuter 
Program 2 assumes first that the surface is a squared sur- 
face and following the same procedure described in Computer 
Program 1, calculates the coupling between dipoles of the 
squared surface. The program has been designed in this way 
in order to be able to apply symmetry. 

DO LOOPS 501 through 507 assign the corresponding 
coupling from the squared surface to ee sehaned surface. 
Then, the same procedure as in Computer Program 1 is applied 
to calculate the right hand side vector, to invert the co- 
efficient matrix, to multiply the inverted matrix by the 
right hand side vector and finally to calculate the current 
density at each dipole. 

The magnitude of the induced current density 1s 


plotted using the three-dimensional plot sub-program PLT3D1 


and is shown in Figure 14. 
38 
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Figure 12. Square Plate (1A by 1A) Magnitude of Induced 
Cibiahareve aise 
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E 
Figure 13. 'L'-Shaped Structure Segmented in 60 Overlapping 
dipoles. Plane Wave Incident at Broadside. 
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'L'-Shaped Structure (1A by 1A) Magnitude of 
Induced Current. 


Figure 14, 
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4. The ‘'T'-Shaped Structure 

The 'T'-shaped structure was segmented in 64 dipoles 
as shown in Figure 15. 

Computer Program 3 listed in Appendix A was used to 
solve the current density distribution for the 'T't-shaped 
structure. The procedure is exactly the same that the one 
used in the 'L'-shaped structure. The differences are in 
the number of dipoles used and the width of the dipole 
me s2 Xk instead of 1/30 A used in the 'L'). 

The current density distribution is shown in Figure 
16 and was plotted using the three-dimensional plot sub- 
program PL3TD1. 

5. The Two-Dimensional Cross 

The two-dimensional cross was segmented in 80 over- 
lapping dipoles as shown in Figure 17. Computer Program 4 
listed in Appendix A was used to solve for its current den- 
Sity distribution. The procedure is exactly the same that 
the one used for the 'L'-shaped structure. The differences 
femeaeche number of dipoles used and the size of the di- 
poles. For the two-dimensional cross D = 1/8 X and w = 
ie 6rd. 

The current density distribution is shown in Figure 
18 and was plotted using the three-dimensional plot sub- 


program PL3TD1. 


Bi. GENERAL CURRENT MODEL 
lier eEnmomecollution of Current distribution in the several 


two-dimensional structures in which two directions of current 
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Figure 15. 'T'-Shaped Structure segmented in 64 Overlapping 
Dipoles. Plane Wave Incident at Broadside. 
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'T'-Shaped Structure (1A by 1A) Magnitude of 
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Induced Current. 


Figure 16. 








VE 
Figure 17. Two-Dimensional Cross Segmented in 80 Overlapping 
7 Dipoles. Plane Wave Incident at Broadside. 
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Two-Dimensional Cross Structure Magnitude 


Of ndweca Current. 


Figure 18. 
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was allowed, the Computer Program written by K. Khemayodhin 
was used. This thesis was concerned only with the design 

of the subroutines to generate the geometry of the structures. 
For all the solutions an electric field of 1201 volts/meter 
incident at broadside was assumed. 

Computer Program 5 as listed in Appendix A was used to 
generate the geometry of all the structures. This program 
generates a square plate of dimensions 1 A by 1 A. The 
difference in the particular structures that has been solved 
1s the current path assumed. Thus, this program can be also 
used for the Unidirectional Current Model by only specifying 
currents in one direction. The reason for not using it in 
the Unidirectional Current Model is the much longer computer 
running time it takes, since symmetry considerations are not 
fiey applied. 

All the results have been plotted using the three~dimen- 
sional plot subprogram PL3TD1. 

ime ithe Rectangular Plate 

The structures of dimensions 1 A by 0.5 A was divided, 
as shown in Figure 19, in 32 surface (monopoles). Eight cur- 
rent paths were assumed in the z-direction and four in the 
y-direction. Figure 20 shows the magnitude of the induced 
current in the direction of the field. FPigure 21 shows the 
magnitude of the cross-polarized induced current. 

ae sauare Plate 

The structure of dimensions 1 A by i A was divided 


in 64 monopoles, as shown in Figure 22. Eight current paths 
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were assumed in the z-direction and eight in the y-direction. 
Figure 23 shows the magnitude of the induced current in the 
direction of the field. Figure 24 shows the efeneeeae of 
the cross-polarized induced current. 
eee Loew Uy =shieaped Structure 

The 'L'-shaped structure was divided in 48 monopoles 
as shown in Figure 25. The dimensions of the Womiventell arm 
tae Dy 0.5 4, and the vertical arm is 0.5 A by 0.5 i. 
Elght current paths were assumed in the z-direction and 
eight in the y-direction. Figure 26 shows the magnitude of 
the induced current in the direction of the field. Figure 
27 shows the magnitude of the cross-polarized induced 


current. 


The 'T'-shaped structure was divided in 48 monopoles 
as shown in Figure 28. The horizontal arm is 1 A by 0.5 4 
and the vertical arm is 0.5 i es eo). eee currant paths 
were assumed in the z-direction and eight in the y-direction. 
Figure 29 shows the magnitude of the induced current in the 
direction of the field. Figure 30 shows the magnitude of 
the cross-polarized induced current. 

5. The Two-Dimensional Cross 

The two-dimensional cross was divided in 48 monopoles 
as shown in Figure 31. The horizontal arm is 1 A by 0.5 A 
and the two vertical arms are 0.5 A by 0.25 A. Eight current 
path were assumed in the z-direction and eight in the y-di- 


mection. Figure 32 shows the magnitude of the current induced 
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Figure 26. 
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magure 27. ‘'L*'-Shaped Structure (1A by 1A) Cross-Polarized 
Current. 
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Figure 32. Two-Dimensional Cross (1A by 1A) Current in 
Ene bareerton. of the Field. 
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in the direction of the field. Figure 33 shows the magnitude 


of the cross-polarized induced current. 
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ieee NUMER ECAR KESULTS 


The three-dimensional plots shown in the preceding section 
give a very good representation of the current density dis- 
aeacion over the two-dimensional structure. However, due 
to the fact that some of the points plotted come from inter- 
polation and that the sub-program PL3TD1 is designed to plot 
only rectangular or squared surfaces, some accuracy is lost. 

In the tables shown in this section, the exact values at 
the specific points are given. 

For the Unidirectional Current Model the value of cur- 
rent 1s at the center of the dipole. For the General Cur- 
rent Method, the value of current is at the edges of the 


monopoles. 
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Table III. Rectangular Plate (Unidirectional Current Model) 
Magnitude of the Induced Current at the Center 
CO Saheewipele. See Figure 9, 
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Table IV. Square Plate (Unidirectional Current Model). 
Waste EMce “Otethestnawced Current at the Center 
of the Dipole. See Figure ll. 
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Table V. 'L'-Shaped Structure (Unidirectional Circuit Model). 


Magnitude of the Induced Current at the Center of 
Pile lpe Mew see Fipure 13. 
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'T'-Shaped Structure (Unidirectional Current 
Model). Magnitude of the Induced Current at 
pew eentereo: tne Dipole, See Figure 15. 
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VV. CONCLUSIONS AND RECOMMENDATIONS 


The results shown for the me cconeules plate (Unidirec- 
tional Current Model) has been compared with those presented 
by Richmond and Wang [Ref. 3]. The results shown for the 
square plate (General Current Model) has been compared with 
those presented by Rahmat-Samii and Mittra [LRef. 4]. Both 
comparisons are seen to be in good agreement. 

A comparison of the results for the Unidirectional Cur- 
rent Model and the General Current Model for similar geo- 
metries, suggest that a solution for the current distribution 
in two-dimensional structures is not complete if the cross- 
polarized current is neglected. This problem is more re- 
markable where the junction problem is present, such as in 
the 'L', 'T' and cross shaped structures. In these struc- 
tures the magnitude of the cross-polarized current reaches 
values of the same order of magnitude as that for the cur- 
rent in the direction of the field, and in some cases as 
for the 'T'-shaped structure and the two-dimensional cross 
the maximum values found for the Brose ens Zed current 
are somewhat larger than those of the current in the direc- 
imeem Of the field. 

The computer programs used in the General Current Model 
G@eenot fully apply symmetry in the calculation of the 
coupling between dipoles. Since this calculation is the 
predominant in computer time consumption, some improvements 


should be investigated in order to be able to solve more 
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complicated geometries without the restriction of "computer 
time." 

A program of surface current measurement for flat-plate 
structures with corner notches, similar to the 'T', 'L' and 
cross structures, should be valuable for validating these 


and other computer codes for two-dimensional surfaces. 
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